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O Rare Kaon decays provide sensitive probes of short distance physics. 

r^i Among the most promising modes are the — t- tt+z/z/ and Kl — )■ n^uu 

^ modes due to their theoretical cleanness and the dedicated experiments, 

^ NA62 and KOTO. The key to the success of New Physics searches in these 

modes is the precision in the standard model prediction. We review the 
^ status and the recent progress of their standard model prediction, and also 

discuss the predictions and New Physics sensitivities of the Kl — )■ ir^i^i^ 
^ and Kl fi~^fi~ decays. 
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1 Introduction 



Within the standard model (SM), short distance (SD) physics accounts for approx- 
imately 90%, 40% and 40% of the branching fractions of the rare Kaon decays 
K — )■ Tipu, Kl — )• 7i^i~^i~ and Kl — )• respectively. These decays are all 

governed by a flavour changing neutral current (FCNC) transition and therefore loop 
induced within the SM. Loop induced processes are special in the sense that heavy 
new particles may contribute to the decay amplitudes on the same footing as SM par- 
ticles. Thus, they represent ideal probes for searching and disentangling New Physics 
(NP). To succeed in isolating NP contributions it is essential to have an accurate 
SM prediction with good theoretical control over both short- and long-distance (LD) 
effects. 

In the following we briefly discuss the general structure of rare Kaon decays and 
to which extend each decay channel can reveal information about degrees of freedom 
beyond the SM. We review the SM prediction of the above-mentioned modes and also 
present the recent electroweak (EW) calculation on the top-quark contribution to the 
two K — i- TTuu decays. 

2 Structure of Rare Kaon Decays 

Rare Kaon decays proceed within the SM via the quark-level transition s — )■ d 
through the box and 7, Z-penguin diagrams of Fig. 1. However, at the Kaons' mass 
scale of approximately 500 MeV the physics is best described by an appropriate 
effective field theory (EFT) in which heavy particles like top-quark, W- and Z-bosons 
are no longer dynamical degrees of freedom. A master amplitude then applies to all 
meson decays [1]: 



(Qi) are matrix elements of effective operators, that cannot be calculated pertur- 
batively, V^km comprise products of CKM matrix elements and Fi are the so-called 
Wilson coefficients, effective coupling constants that parametrise effects of SD physics. 
The sum in Eq. (1) extends over all possible operators Qi generated by the SM or 
a given NP model and similarly over contributions from different up-type quarks in 
the loops of Fig. 1. The EFT approach allows to separate SD from LD physics and 
therefore to treat these contributions independently. While the matrix elements of 
semileptonic decays can be extracted from experimental data, the Wilson coefficients 
depend on the high-energy model and can be calculated perturbatively within the 



Two points should be stressed regarding Eq. (1). First, are process independent 
universal loop functions within the SM, meaning that they are the same for B-, D- 
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SM. 
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Figure 1: 7, Z penguins and W boxes contributing to the s — >• d transition at the one- loop 
level. In the case I = v no 7-penguin contributes. 



and K-meson decays. The resulting correlation of observables in the SM tests the 
universality of the loop functions in a model independent way (see talk by D. Straub 
[2]). Secondly, NP can not only change the Wilson coefficients, but can also generate 
additional operators contributing to the decay amplitude. This would signal a depar- 
ture from Minimal Flavour Violation [3]. Both NP effects can affect the branching 
ratios of rare Kaon decays. 

The master amplitude of Eq. (1) also incorporates the relative magnitude of each 
contribution. Its size depends on both the CKM elements, conveniently parametrised 
in terms of A = \Vus\ ~ 0.22, and on the loop functions Fj's. Within the SM, 
low-energy contributions are suppressed with respect to the top contribution due to 
the Glashow-Iliopoulos-Maiani (GIM) mechanism [4]. When SD photon penguins 
contribute, this suppression is only logarithmic, but in their absence (e.g. in the 

TTvu decays) the loop functions are proportional to ml^p, "^charm' ■'^qcd 
the up-quark, and a hard quadratic GIM is at play. 

In addition, the CKM factor of the top-quark contribution At = V^*Vtd is propor- 
tional to A^ and therefore highly suppressed. This strong suppression, absent in B 
and D decays, can imply a significant charm-contribution depending on the mode, 
but also renders rare Kaon decays sensitive to even very small deviations from the 
CKM picture of CP violation. 



3 K'^ TT^uu and Kl TT^ff 

The weak SM effective Hamiltonian for the two K — )■ vrz/z/ decays reads [5] 

^eff= ^^^ ."2^ V {\X + \Xt){-SLl,dL){uiLrUlL) + h.C. (2) 

V2 27r sm 6w , 

and involves to a good approximation only one effective operator below fi = fic ■ 

The LD part of both decays is known with high precision. The matrix elements are 
extracted from K^^ decays including isospin-breaking (k+'^) and long-distance QED 
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Figure 2: Left: comparison of the top-quark Wilson coefficient in the MS and on-shell 
schemes at LO, NLO large-rrif and NLO. The grey sohd band, hashed band, and sohd red 
band represent the estimated uncertainties of LO, NLO large-m^, and NLO result as a 
function of the Higgs mass, respectively. Right: the uncertainty due to EW renormalisation 
scale at LO and NLO. Again the grey and red band represent the estimated uncertainty 
before and after the NLO calculation, respectively. 



effects [6], whereas also effects from light quarks and higher dimensional operators 
{SPc,u) have been estimated using chiral perturbation theory (ChPT) [7]. 

Regarding the small uncertainties from LD effects, higher order corrections on 
the SD parts become essential for a precise SM prediction. The SD charm contri- 
bution, denoted by Pc = il-^'^ + l-^^)^ negligible for — )■ vr^z/z/, but not for 
— )■ n'^uij, where it amounts to approximately 30% of its branching ratio. So far 
Pc is known up to next-to-next-to-leading order (NNLO) in QCD [8] and also the 
next-to-leading order (NLO) EW corrections have been calculated [9]. NLO QCD 
corrections [5, 10] are known for the top-quark contribution, but until recently the 
NLO EW corrections were only known in the large-m^ limit [11]. However, these were 
known to poorly approximate the full EW corrections (see Fig. 2) [11]. Therefore, the 
renormalisation scheme of the EW input parameters a, sinOw, Myy, Mj, appearing at 
LO in Eq. (2), was not clear and accounted for an approximately 2% uncertainty in 
Xt, which scales quadratically in the uncertainty of the branching ratios. 

Our recent two-loop calculation of the full NLO EW corrections on Xt [12] sub- 
stantially reduced the renormalisation scheme uncertainty to approximately 0.2% and 
also removed the remaining scale dependence. The comparison of LO and NLO re- 
sults in the MS and on-shell scheme is illustrated in Fig. 2 together with the remaining 
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Figure 3: 7- and 77- penguins contributing to the LD part of the Kl — )• ir and 
Kl — ^ M^A*^ branching ratios. 

scale dependence. 

The theoretical predictions and the current experimental values for the branching 
ratios then read: 

= (7.81^°:?? ± 0.29) X 10-1^12] BR-?^^,^, = (1.73^1;^^) x 10-''[13] 
BR*£°_.o., = (2.43^°:^? ± 0.06) X 10-" [12] BR-^^,,^, < 6.7 x 10"^ [14] 

The first error in the theory prediction is parametric, while the second summarises the 
remaining theory uncertainty. For the charged mode the detailed main parametric 
contributions are {Vcb '■ 56%, p : 21%, nic : 8%, rrit : 6%, f] : 4%, : 3%, sin^ 9]v : 1%), 
while the main theoretic contributions are (5Pc,« : 46%, Xt(QCD) : 24%, Pc : 20%, 
: 7%, Xt(EW) : 3%), respectively. Similarly for the neutral mode, the parametric 
uncertainties are {Vcb '■ 54%, fj : 39%, rrit : 6%) and the theoretical are (Xi(QCD) : 
73%, : 18%, Xi(EW) : 8%, SP^^u : 1%), respectively. 

In the two Kl — > 7r°£"^£~ modes three different contributions compete in magnitude: 
direct- CP- violating (DCPV), indirect-CP- violating (ICPV) and CP-conserving (CPC) 
contribution. Long distance effects are much larger here. Still, the importance of 
these decay channels lies in their sensitivity to SD contributions from more than one 
operator. Let us now look at the three distinct contributions. 

DCPV: high-energy top and charm contributions generate both vector, Qjv = 
{sd)v{i£)v, and axial-vector, Qta = {sd)v{i^)A, dimension-six operators. Their Wil- 
son coefficients are known to NLO QCD [3] and their matrix elements are also ex- 
tracted with a few per mil precision from Ki^ decays [6] . One difference between the 
electronic and the muonic channels is the suppression of the electronic axial-vector 
matrix element due to the smallness of the electron mass. 

ICPV: in this case the Kl decays through its small CP even component via the 
7-penguin in Fig. 3. The size of this smaller component is described by the parameter 
ex- We can therefore relate this contribution to the decay Ks{^ Ki) — )■ 7r°£+£~ using 
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ChPT. The light meson loops are subleading and the amplitude is dominated by the 
chiral counterterm as [15], whose 20% uncertainty completely dominates the error for 
the Ki — i- 7i^i^i~ rates [16]. Although the absolute value of as has been extracted 
from NA48 measurements, its sign remains unknown. Therefore, it is still an open 
question whether the ICPV contribution interferes constructively or destructively 
with the DCPV contribution. Apart from a more precise Ks — )■ 7i^i~^i~ measurement, 
the sign of as can also be determined from the measurement of the integrated forward- 
backward asymmetry Ap^ in — 7r°yU~^/i~ [16]. 

CPC: is a purely LD effect described by the 77-penguin ChPT diagram in Fig. 3. 
The contribution is helicity suppressed for the electronic mode, a further difference 
between the two modes. For the muonic mode the meson loops are finite at O(p^), 
while higher order effects of 0{p^) have been partially estimated from the Kl — )■ 7r°77 
ratio, reducing the uncertainty to approximately 30% [17]. 

The current SM predictions and experimental values for the branching ratios for 
the two decays then read: 

BRZ%.o,,,- = 1.41;°:^^ (0.95l°|?) X 10-"[16] BR-^^o,+,- < 3.8 x 10-i°[18] 
BR*^^r-..Oe+e- = 3.54l°|« (1.56^^9) X 10-^^16] BR^^^o,.,- < 2.8 x lO'^" [19] 

Theory predictions correspond to constructive interference of DCPV with ICPV con- 
tribution, while predictions in brackets to destructive interference. 

5 Kl ^ 

Kl — ^ /U'^/i" is a CP conserving decay, whose SD sensitivity is obscured by large LD 
effects. The top contribution is known up to NLO QCD [20] and the charm contri- 
bution up to NNLO QCD [21]. Indirect CP violation is negligible in this channel. 
However, in contrast to the Kl — t- 7r°yU+/i~ case, the LD 77-penguin is not finite 
(Fig. 3). The absorptive part of the 77-loop, extracted from Kl — )■ 77, accounts to 
approximately 98% of the measured branching ratio, BR^^_^^+^_ = 6.84(11) x 10~^ 
[22] , while dispersive part of the 77-contribution is divergent in ChPT and can there- 
fore be estimated from Kl — )■ 7*7* only with a large uncertainty [23]. On top of these 
difficulties, the 77-contribution interferes with an unknown sign with the SD contri- 
bution, which further increases the theoretical uncertainty of the channel. Therefore, 
the interesting SD part of the decay mode cannot be extracted accurately. Neverthe- 
less, useful NP constraints can be derived using its measured branching ratio. 

6 New Physics and Outlook 

The theoretical cleanness of the K -nvv modes promotes both to excellent probes 
of NP, especially in view of the dedicated experiments, NA62 at CERN and KOTO at 
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JPARC, which aim at measuring the charged and neutral mode, respectively, with an 
expected accuracy of 15%. These decays have therefore been studied extensively in 
models beyond the SM. A summary of predictions in different models is presented in 
[24]. This analysis points out the possibility of large effects in both decay modes and 
also illustrates the correlation between the predictions of the two decays in a large 
class of NP models. 

NP effects on the branching ratios of the — Tr^f"*"^" modes have also been 
considered, in a model-independent way. A measurement of the decays could test con- 
tributions from operators not generated by the SM. Effects from scalar, pseudoscalar, 
tensor, and pseudotensor operators, generated by NP with or without helicity sup- 
pression, have been studied and can be disentangled by a more precise measurement 
of both modes [16]. 

To conclude, rare Kaon decays are very clean and sensitive probes of NP. The 
observation of deviations from the SM by the dedicated experiments would be a clear 
NP signal. Also, the use of both predictions and measurements as constraints of NP 
can shed light on NP patterns and the structure of flavour violation. 
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